Formal Valence, d Occupation, and Charge-Order Transitions 
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While the formal valence and charge state concepts have been tremendously important in materials 
physics and chemistry, their very loose connection to actual charge leads to uncertainties in modeling 
behavior and interpreting data. We point out, taking several transition metal oxides (La2VCu06, 
YMO3, CaFeC>3, AgNi02, V4O7) as examples, that while dividing the crystal charge into atomic 
contributions is an ill-posed activity, the 3d occupation of a cation (and more particularly, differences) 
is readily available in first principles calculations. We discuss these examples, which include distinct 
charge states and charge-order (or disproportionation) systems, where different "charge states" of 
cations have identical 3d orbital occupation. Implications for theoretical modeling of such charge 
states and charge-ordering mechanisms are discussed. 



Spin ordering, and often orbital ordering, is nor- 
mally unambiguous, as these properties are subject 
to direct observation by magnetic and spectroscopic 
measurements, respectively. Charge ordering (CO) 
and the actual charge of an ion is rarely measured 
directly, and the formal charge of an ion in the solid 
state can be a point of confusion and contention. 
Valence, oxidation number, and formal charge are 
concepts borrowed from chemistry, where it is em- 
phasized they do not represent actual chargefH Q 
and have even been labeled hypothetical. [l[ As the 
interplay between spin, charge, orbital, and lattice 
degrees of freedom become more closely watched Q 
and acknowledged to be a complex phenomenon, dis- 
proportionation and CO have become entrenched 
as the explanation of several high profile metal- 
insulator transitions (MIT). The possibility that CO 
in the charge transfer regime is associated with the 
oxygen sublatticc, with negligible participation of 
the metal, has been raised[J| and considered as an 
alternative. Q 

Charge density is a physical observable of con- 
densed matter, and the desire to assign charge to 
atoms has evident pedagogical value, so theoretical 
approaches have been devised to share it amongst 
constituent nuclei. Mulliken charge population, 
which socializes shared charge (divides it evenly be- 
tween overlapping orbitals) is notoriously sensitive 
to the local orbital basis set that is required to spec- 
ify it. Born effective charges are dynamical proper- 
ties and are often quite different from any conceiv- 
able formal charge or actual charge. Integrations 
over various volumes have been used a great deal, 
but dividing the static crystal charge density into 
atomic contributions is, undeniably, an ill-defined 
activity. 

A possibility that has not been utilized is that, 
taking 3d oxides as an example, there is a directly 
relevant metric that is well defined: the d occupa- 
tion rid- This quantity is in fact what the physical 
picture of formal charge or oxidation state brings to 
mind. 3d cations, in their various environments and 



charge states, have maxima in their spherically aver- 
aged radial density p(r) in the range 0.6-0.9 a . At 
this short distance from the nucleus, the only other 
contribution to the density is the core contribution, 
which can be subtracted out and is unchanged dur- 
ing chemical processes or CO. Most relevant to the 
understanding of CO-driven transitions and dispro- 
portionation is the (actual or relative) difference in 
3d occupations Arid, which is given directly, without 
any integration, by the difference in the radial 3d 
densities at their peaks, where there are no compet- 
ing orbital occupations to confuse charge counting. 
This specifically defined 3d occupation differences 
provides a basis for building a faithful picture of 
CO and of characterizing formal valence differences 
more realistically. We consider our computational 
results^, 0] for a selection of systems, then discuss 
some of the implications. 

La2VCuQ 6 (LVCO) is a double perovskite com- 
pound providing a vivid and illustrative example. 
Our earlier studyQ revealed two competing config- 
urations for the ground state. Using conventional 
identifications, one is the V 4+ d 1 , Cu 2+ d 9 magnetic 
configuration (with bands shown in Fig. [T|) identi- 
fied as such because (1) there is one band of strong 

V d character occupied and one band of strong Cu d 
character unoccupied, and (2) the moments on both 

V and Cu, 0.7 [Ib, are representative of many cases 
of spin-half moments reduced by hybridization with 
O 2p orbitals. The other configuration is the non- 
magnetic d° — d 10 band insulator: all Cu d bands 
are occupied, all V d bands are unoccupied - a con- 
ventional ionic band insulator in all respects. The 
identification of formal valence (or oxidation state) 
is crystal clear. 

The radial charge densities of V and of Cu for both 
configurations reveal an unsettling feature: the ac- 
tual 3d occupations nd of each of these V and Cu 
ions are identical for both configurations, in spite of 
the unit difference in their formal charges. (Identical 
in this paper means to better than 0.5% ~ O.Ole - , 
in terms of the differences of charge density at their 
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FIG. 1: (Color online) Top: bands near the Fermi en- 
ergy/bandgap in the d — d 9 magnetic, nearly Mott in- 
sulating, configuration of La2VCu06- The d xy -up band 
is correlation-split off from the other two ti g bands and 
fully occupied. Bottom: the Cu fatbands for the same 
system, showing one unoccupied Cu minority d x 2_ y 2 
band correlation-split from the d z i band. The other d 
bands fall outside this energy range. 



peaks.) Thus ions with no real difference in 3d occu- 
pation can behave as if they comprise charge states 
differing by unity. Changes in spin-orbital occupa- 
tions, which quantify spin, orbital, and charge differ- 
ences between the two states, can be quantified by 
the LDA+U spin-orbital occupations. For the V d 1 
d xy (Jahn- Teller split) orbital, the majority-minority 
difference is 0.70, which accounts for all of the mo- 
ment. The difference of 0.65 between d xy and each 
of the other ti g characterizes the Jahn- Teller dis- 
tortion. The increase in charge of the d xy orbitals 
(both spins), 0.55, compared to the d° state, is ab- 
sorbed more or less uniformly from all other (nom- 
inally unoccupied) spin orbitals. Similarly for Cu, 
the d 9 hole results from a difference of charge in the 
minority d x 2_ y 2 orbital of 0.6, with the other hole 
charge being distributed nearly uniformly over the 
other nine (nominally but not actually fully occu- 
pied) spin-orbitals. In both cases the moment arises 
entirely from the single magnetic orbital as the sim- 
ple picture would suggest, while all other orbitals 
are unpolarized. This happens, conspicuously, with 
no change in nd for either V or Cu. Charge is redis- 



tributed to one orbital from the others, and strongly 
spin-imbalanced within that orbital. Even with in- 
sulators with "obvious" charge states, 3d orbital oc- 
cupations can range over the values [0,1]. 

We look at additional cases before addressing 
some of the implications. 

Rare earth (TV) nickelates 7?Ni03 display a first 
order structural and MIT of great current interest. 
The Pbnm (GdFeOs structure) — > P2\/n transfor- 
mation results in a large NilOe and a small Ni20g 
octahedron, with Ni-0 distances of 2.015±0.015 
A and 1.915±0.025 A, respectively, that are not 
otherwise strongly distorted; see the inset of Fig. [2] 
At a temperature that varies smoothly from 600K to 
300K with increasing TZ ionic radius, the resistivity 
of these nickelates drops sharply. 0, [l(| We focus on 
YNIO3; with its small ionic radius, it is one of the 
more strongly distorted members, and the resulting 
narrowed bandwidths make it more prone to strong 
correlation and CO tendencies. [H[ Structural changes 
at the MIT have been studied extensively,]^, fill — 



14| which together with x-ray absorption spectral 



splittings jl5Hl7j ] have been interpreted in terms of 
charge disproportionation (or CO) 2Ni 3+ -> Ni 3+<5 
+ m 3 ~ s , with 5 w 0.3 for YNi0 3 .[l5| 

This MIT in the nickelates has been recognized as 
paradigmatic by theorists. Mizokawa et al. mod- 
eled this system 18] with a multiband Hartree-Fock 
model in the charge-transfer regime and found evi- 
dence for CO on the oxygen sublattice for larger TZ 
cations, but concluded that YMO3 was representa- 
tive of a CO transition on the Ni sites. Mazin et 
al. [H surveyed the competition between Jahn- Teller 
distortion of the d 7 ion and CO and also concluded 
that YMO3 is a prime example of a CO d 6 + d s 
system. Lee et al. have investigatedflij a two band 
model for this system with a CO interaction in mean 
field, emphasizing CO effects. On the other hand, 
Yamamoto and Fujiwara[2(| reported a very small 
(^0.03 e~) density functional based charge differ- 
ence. 

For the assumed (for simplicity) ferromagnetic or- 
der the calculated Nil and Ni2 moments arc 1.4 
and 0.65 fiB respectively for YN1O3 and several 
other members of this class, so these values are 
not sensitive to the magnitude of the distortion. 
They coincide with the values obtained from neutron 
diffraction, [lH 1.4(1) and 0.7(1) /is respectively, in 
the magnetically ordered phase. It is intriguing that 
the same moments were obtained in fully relaxed 
LaNiOs/LaAlOs monolayer superlattices.[21| 

The 3d occupations, obtained as above directly 
from the maximum in the radial charge density plots 
in Fig. [21 arc identical for Nil, Ni2, and the single 
Ni site in the high temperature phase: there is no 
3c? charge transfer, or disproportionation, across the 
transition. The majority and minority radial densi- 
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FIG. 2: (Color online) Radial charge density (upper 
curve) of YNi0 3 for Pbnm Ni and P2 1 /n Nil and Ni2, 
showing there is no difference at the peak, which reflects 
the 3d occupation of the ion; a small difference shows 
up near the sphere boundary. The spin decompositions 
give easily visible differences. The vertical lines at the 
bottom right indicate conventional Ni + , Ni 3+ , and Ni 2+ 
ionic radii, which have no relation to the (unvarying) 3d 
occupation. Inset: Structure of the broken symmetry 
P2\jn phase, showing the rotation in the a — b plane 
and tilting along the c axis of the NiC>6 octahedra (Ni 
is inside) and the (ir, ir, ir) ordering of the Nil and Ni2 
octahedra. 



for in any simple model. In spite of identical 3d 
charges, the Nil and Ni2 core energies differ by up 
to 1.5 eV. 

CaFeC>3, another perovskite that displays the same 
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FIG. 3: (Color online) Spin-decomposed Ni ti g and e g 
density of states for the Nil and Ni2 ions in the insulat- 
ing P2\/n broken symmetry phase. The hashed regions 
illustrate the spectral origin of the enhanced moment of 
Nil relative to Ni2. The horizontal arrows illustrate the 
large difference in spin splittings, the result of the com- 
bination of Hund's coupling and Coulomb U= 5.7 eV. 



ties and integrated charges of course differ (see Fig. 
^ as they must to give the moment, but the total 
3d occupation is inflexible. This constancy of the 
3d occupation across the transition, and equality for 
Nil and Ni2, is inconsistent with microscopic dispro- 
portionation. 

To illustrate the spin-orbital spectral density re- 
distribution, the projected densities of states are 
shown in Fig. [3] All ti g states are filled and irrel- 
evant. The e g spectral distribution is non-intuitive: 
weight from -5 eV spin-down is transferred to -1 eV 
spin-up. The majority e g states just below the gap 
are strongly Nil in character, while the unoccupied 
bands just above the gap are primarily Ni2. Such 
behavior is expected for different charge states, sim- 
ilarly to the behavior in LVCO above; however, the 
total 3d occupation is identical. 

The main differences between Nil and Ni2 show 
up in the unoccupied e g states: the Nil spin split- 
ting is 3.5 eV, a reflection of the on-site repulsion 
that opens the Mott gap in the majority e g states, 
rather than Hund's exchange splitting. The origin of 
the Ni2 moment is murky, not identifiable with any 
occupied spectral density peak. Note that in a Ni 2+ 
+ Ni 4+ CO picture, Ni2 would be nonmagnetic. Not 
only is this calculated behavior not consistent with a 
CO picture, it involves redistribution not accounted 



Pbnm — > P2\/n structural change at Tmi as the 
nickelates, is also cxplained[22[ in CO language that 
invokes the unusually high (penta)valent state Fe 5+ . 
Analogously to YMO3, we obtain identical 3d oc- 
cupations for Fel and Fe2 ions. Quantum chemi- 
cal embedded cluster calculations [23j and LDA+U 
studies [UEll had noted that the Fe charge in both 
"disproportionatcd" sites differed little, but neither 
quantified the occupation as we have for YMO3 and 
CaFeOs. The pentavalent state of Fe has most often 
been identified from Mossbauer isomer shift data, 
but Sadoc et a/.Q concluded the difference in iso- 
mer shift is primarily a measure of the covalency 
(Fe-0 distance) rather than any real charge on Fe. 

AgNi02, a triangular, magnetically frustrated lat- 
tice compound with nominal Ni 3+ ions, undergoes 
a structural transition at 365 K although remain- 
ing metallic. 27- 3^1 Three inequivalent Ni sites arise, 
with a high spin Nil ion in an enlarged octahedron 
and two low spin Ni2, Ni3 = Ni2,3 ions in small 
octahedra. Based on the structural changes (which 
were quantified in terms of bond valence sums), the 
magnetic moments, and resonant x-ray scattering 
that confirms a calculated ~1 eV difference in core 
level energies between Nil and Ni2,3, this transi- 
tion has been welcomed as the first realization of 
such a highly unusual 3e g — > e 2 g + 2eg 5 type of CO. 
Furthermore, using the charge difference per unit 



core level splitting of 0.66 e/eV led to an inferred 
charge disproportionation of ^1.65e, i.e. Nil 2+ + 2 
Ni2,3 3 5+ . We have reproduced several of the first 
principles results [271 |30T ] that were used to support 
CO. The calculations give a large moment (> 
on high-spin Nil and very weak moments (~0.1 /is) 
on low-spin Ni2,3 ions. We find, as in the cases 
above, that rid for the three sites are identical. More- 
over, our calculated core level differences, 0.6-0.8 eV, 
are roughly consistent with reported values [3Cj| (~1 
eV). 

V4O7 represents another oxide currently ex- 
plained by a CO-driven MIT. It is structurally more 
involved, but first principles calculations of mo- 
ments and geometries again have produced several 
results corroborating the experimental data[3l|, [32| 
and were used to support CO into V 3+ and V 4+ 
charge states on specific sites. As in the instances 
above, we find no differences in rid', the occupations 
are indistinguishable. The site energy differences, 
measured by differences in Is, 2s, 2p core levels, dif- 
fer by 0.9-1.2 eV for two sites, similar to the nickc- 
lates. The interplay of orbital order, structural dis- 
tortions, and possible spin-singlet formation of half 
of the V ions provide a rich array of degrees of free- 
dom, which can operate without need for dispropor- 
tionation. 

Implications. We have established that, for sev- 
eral instances of CO transition systems as well as 
for the two self-evident charge states of LVCO, there 
is no difference in the 3c? occupations for the dif- 
ferent "charge states" that have been used to cat- 
egorize their behavior. Such identification is possi- 
ble because a choice of a region for integration is 
avoided; the peak charge region rather than tails of 
orbitals are used in the identification. This find- 
ing of constancy sharpens several reports of "small 
charge differences" between differing charge states 
(viz. Luo et al. 33[ for doped manganitcs; Haldanc 
and Anderson [341] in a multi-orbital Anderson model, 
and Raebiger et aZ.[35[ from DFT calculations for 
TM impurities in semiconductors; Yamamoto and 
Fuiiwaraj2oT| and also Park et al. 36] for nickelates). 

We see two primary implications: (1) the concep- 
tual basis underlying a substantial aspect of transi- 
tion metal physics is misleading, and (2) modeling 
of structural and electronic transitions has, at least 
in several conspicuous cases, incorporated the wrong 
mechanisms by invoking inactive degrees of freedom. 
Actual cases of CO very likely do exist, but the bur- 
den of proof has shifted. 

For these CO systems, the constancy of rid sug- 
gests that Ud is too large to allow change in occupa- 
tion rid in or near the ground state (in the cases we 
discuss, and similar ones). In insulators the charge 
is more physically pictured in terms of (fully occu- 
pied) Wannier functions (WFs) than in terms of am- 



biguous populations of atomic orbitals, making them 
appear to be inviting. However, WFs are far from 
unique and, like molecular orbitals, WFs contain 
charge that cannot objectively be assigned to one 
atom or another, so a WF viewpoint is not promis- 
ing. 

A broader implication is that modeling of cou- 
pled structural and electronic transitions in terms 
of charges [H [l9| from atomic- like orbitals must be 
treated with caution: charge fluctuations in these 
systems are too high in energy to comprise a relevant 
degree of freedom. The important energy differences 
are characterized in terms of differences in hopping 
amplitudes, anion-cation distances, and (not recog- 
nized in most models) resulting changes in site ener- 
gies, as well as very important Hund's rule energies. 
Models that try to parametrize (for example) Nil- 
Ni2 differences by on-site charge will not be treat- 
ing the relevant microscopic degrees of freedom. CO 
on the oxygen sublatticejj, [B| may also be problcm- 
matic. 

Charge states of ions serve to specify the occu- 
pations of spin-orbitals. The essential degrees of 
freedom in determining this popular characteriza- 
tion, which professes to be quantitative, are the 
spin-orbital occupations, not as determined from the 
(real) density matrix but rather from the site sym- 
metry, crystal symmetry, and the local moment. The 
LVCO example illustrates vividly how two different 
charge states, for both highly charged V and moder- 
ately charged Cu, can be represented by integer oc- 
cupation of different numbers of orbitals while there 
is no change in rid- "Charge state" projects onto in- 
tegrally occupied orbitals, while the distribution of 
real charge is strongly non-integral and often non- 
intuitive. These projections are backed up by the 
number of occupied spin-polarized bands (an inte- 
ger), by the (discrete) local symmetry (JT distor- 
tion), by the local moment (with its quantization 
smeared by hybridization), and by the atomic radii, 
but each one of these characterizations is extremely 
flexible with a given amount of 3c? charge. 

More specifically to CO systems, the ionic envi- 
ronment in the high symmetry phase requires closer 
scrutiny. In both the nickelates and in V4O7 there is 
evidence of distinct metal sites above the transition, 
in the (on average) symmetric phase, and the struc- 
tural similarities of CaFeOs to i?M03 suggest simi- 
lar behavior there. For nickelates, x-ray absorption 
spectra [H, 17 1 reveal that local signatures of Nil and 
Ni2 sites persist continuously across the MIT, and 
both sites also remain when driven across the phase 
boundary by pressure. 37J As we have shown, the co- 
ordination alone (i.e. with identical rid) accounts for 
on-site energy differences of ~1 eV in spectra that 
have often been used to support disproportionation. 
The MITs in some of these materials may be pri- 



marily order-disorder type; the onset of long-range 
order in nickelates results in carrier localization and 
gap formation, ergo a MIT but one unrelated to CO. 

We propose therefore that "charge order" should 
be used as the name, hence the interpretation, of a 
phase transition only if an objective, relevant charge 
difference is the likely mechanism; otherwise, the 
underlying mechanisms should be identified. For- 
mal developments may be useful; for example, Jiang 
et al. have provided a specification's] of integer 
charges in an insulator that they propose as oxi- 
dation states (which are identical to charge states 
in metal oxides.) Based on integration over a con- 
figuration space path of the dynamic Born effective 
charge, their expression assigns (in principle) an in- 
teger charge to each atom in any insulator. Notably, 
their specification does not refer to 3d charge ex- 
plicitly and furthermore depends explicitly on dy- 
namical effects (electron response to ion motion). 
Also, many CO interpretations only hold water if 
the supposed charge difference ±5 is much smaller 
than unity (5 ~ 0.3 for the nickelates). More expe- 
rience will be needed to learn how best to interpret 
their definition. 
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